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In continuation of a previous paper a close connection between Feynman propagators and a particular 
Levy stochastic process is established. The approach can be easily applied to the Standard Model 
SUc(3) x SUl(2) x (7(1) providing qualitative interesting results. Quantitative results, compatible 
with experimental data, are obtained in the case of neutrinos. 
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i ,! I. INTRODUCTION 

^. 

In a previous note [l[ we showed that it is possible to derive and modify the relativistic Feynman propagator of 
Q-i a free (forceless) particle (fermion spin i, boson spin and 1) on the basis of Levy stochastic processes Indeed 
we derived the relevant connection between the kinetic energy E$ of a free particle and the logarithmic characteristic 
i ' v( u ) °f a Levy process (see Appendix |A"| : 

E Q = -™ V (u) (1) 



. a being an action, r a time parameter and u = p/a with p the particle momentum. The choice r)(u) = 1 — ^/l + i 
C^l " leads to the well known relativistic formula 



q{ ■ En = cy/m 2 c 2 + \p\ 2 - mc 2 (2) 

o: 

£S| , once the following identifications arc done 



a 2 a 

— = mc , a = — , p 

t mc 



Furthermore it was proved that (f(u) = e^™' is the characteristic function of an infinitely divisible probability 
distribution as it is required for having a Levy process Q. At this point we proposed a modification of equation @ 
that preserves the infinite divisibility and produces an underlying Levy stochastic process Such modification 
consists of inserting a function f(x) with x = p 2 /m 2 c 2 (p 2 = E 2 /c 2 — p 2 , E = E + mc 2 ) into the square root of the 
energy-momentum formula 

E = ™\ll+^ + f(4^) (3) 



The aspects of f(x) are analyzed later in this note; now we observe that the total energy E is defined in the equation ([3]) 
implicitly because of the argument of f(x) containing E (p 2 = E 2 jc 2 — p 2 ). 
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We can interpret, in quantum mechanics, equation (J3j) as an energy equation (Hamiltonian) that acts on a wave 
function ip(x , t) by using the usual prescription E — > ihd t and p — > — £/iV with a = h 



ihdt^p(x,t) = mc Wl 



K 2 V 2 



+ f 



h 2 a 



In the usual way one then achieves the modified Klein-Gordon and Dirac equations (h = c = 1) 

□ - m 2 / f ) - m 2 ^ = 

i\) = 

whereas the corresponding Feynman propagators verify the inhomogeneous equations 

'□ 2 N 



i 7fl 9"- m Wl + /(^ 



□2 - m 2 / 



in/nd% - m* 1 + f 



\m'- 

In momentum space 1 they have the following representation 

£kg{p 2 ) = 

ic d (p 2 ) = 



JCkg(2\1) = S 4 (2\l) 
Kn(2\l) = *<5 4 (2|1) 

1 



p2 _ m 2 ^ _|_ /(j? 2 /m 2 )] + te 
1 



l^Pn ~ wiyl + f(p 2 /m 2 ) + ie 



(4) 

(5) 
(6) 

(7) 

(8) 

(9) 
(10) 



We note that the equations 0, (JSJ), © and (fTU)) coincide with the well known formulae (in absence of interaction) if 
/Or) -> 0. 



II. MODIFIED FEYNMAN'S GRAPHS 



Let us now consider the Feynman pcrturbative contributions in renormalized field theories and focus our attention 
on the self-energy terms in QED, QCD and in SM SUc(3) x SUl(2) x £7(1). The amplitude for a fermion that 
propagates from the vertex X to Y , if expanded, looks as follows (the numerator is not essential for our purposes) 



A = A^+A^+A^ + 



with zero order 



A<°> = Y- 



l^Pn ~ fHv/1 + f(p 2 /m 2 ) + ie 



-X 



(11) 



(12) 



If the fermion emits and reabsorbs a virtual boson with mass M (M = Mz,M w ± for weak interactions, M = for 
photons or gluons) we have 



a^ = Yy 



1 



-c- 



1 



— myl + f(p 2 /m 2 ) + ie r y v p u — myl + f(p 2 /m 2 ) + i 



(13) 



1 The connections between the wave functions and their propagators of the equations J7J and JHJ are then 

V>(2) = J[^(l)d lfl K KG (2\l)-di^(l)K KG (2\l)]N^(l)d 3 Vi 
V>(2) = y?C D (2|l)7^^(l)^(l)d 3 V'i 
Nfj,(l) being the inward normal to the surface at the point 1. 
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where 

C = ( ^9 2 s d A k == = ~ 2 ~ 2 

7 [(p - fc) 2 - M 2 ] [7" fc M - m y/l + f (k 2 /m 2 )} + le 

and g s is the electro- weak renormalized coupling or strong quark gluon coupling Q . 

III. PHENOMENOLOGY 

For M = (gluon renormalized mass, see equation (1141) ) we approximate the exact A(p 2 ) in a simple way by adding 
an infinite number of Feynman graphs that contain only one gluon at any time: 



L> = 7/-mv/l + /(p 2 /m 2 ) (16) 

All other contributions that we neglect here contain two or more virtual gluons simultaneously in the intermediate 
states: they are supposed to be less relevant starting with terms of order g . Hence from equation (|1 5|) we may obtain: 



A(p 2 ) = Y * X = Y J_ X (17) 

7-p-mvl + / -C j.p — my/1 + f - Arf - B 



The search of poles of A(p ) ~ A(p ) leads to the equation 



p 2 (l - A) 2 - m 



B( P 2 ) 



m 



2 



(18) 



Under the approximation A(p 2 ) ~ A(m 2 ) and B{p 2 ) ~ B(m 2 ) that follow from the very reasonable assumption 
A< 1, i? <C myl + f(p 2 /m 2 ) we obtain the equation 



p 2 



1 - i(m 2 ) 



(19) 



which for g s —> 0, A, — > gives back the classical equation 



4=i+/f4) 



Note that equation (f2U| provides poles of p 2 at zero order (no coupling) . More specifically we have 

f(x) ~ A + Aix + A 2 a: 2 + A 3 x 3 (21) 
x-/(x)-l = -\s(x-xo)(x-x+)(x-x-) (22) 

so that from /(l) = we get Ao = — Ai — A 2 — A3 and (for more details see Appendix |B|) 

2 2 

= 1 x + = —| a;_ = —| 23) 
mf mf 

We now remind that the field theories we are considering here (QED, QCD and SUc(3) x SUl(2) x are 
renormalizable and consequently the function f(x) can be looked at as a smooth cut-off that regularizes the per- 
turbative terms. More precisely the integral C becomes finite (integrand convergent) under the assumption that the 
function f(x) is analytic and appears either as a series expansion or as a polynomial of third degree at least. In 
the latter case one obtains three poles in the (zero order) fermionic propagator and obviously three physical masses 
for appropriate values of its coefficients. We may relate this result with the three families of fundamental particles 
(quarks) thus describing the flavour phenomenon: more explicitly two propagators (charge — ^ and +|) for quarks, 
one for charged leptons (e - , t~) and one for neutrinos (y e , u^, v T ). However the numerical values of these masses 
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cannot be considered valid if compared with with experimental data, except for particular cases where the strong and 
e.m. interactions are absent (see neutrinos in the next section). 

We now consider the Feymnam propagator for fundamental bosons (gluons, W , Z. Higgs). Again we develop it in 
perturbation theory and discover that the most significant contribution that must be added to the zero order derives 
from the fcrmion-antifcrmion loops. To be precise one has the zero order propagator that connects a fermion vertex 
7 M with another j v : A = + A^ + ... 



A (o) 



P 



M 2 



9^ + (e - 1) 



P 



£ 2 M 2 



where £ = 1 in t'Hooft-Feynman gauge, £ = oo in unitary gauge and £ = in Landau gauge 



(24) 



A W = 47T3 2 / -^-Tr 
y J (2tt) 4 



q-p 



1 + / 



(<?-p) 2 



--lu- 



1 + / 



(25) 



IV. NEUTRINOS 



Within the scenario of weak interactions, neutrinos need a special attention. Indeed they arc considered to be 
masslcss in the Standard Model but experiments have shown that they have masses. These masses are much smaller 
than those of other fundamental particles. Neutrinos are produced by weak interactions in definite flavor states {y e , 
Vp, v T ) which are not stationary mass (energy) states and one expects mixing and oscillations between the various 
flavor states. Both mixingand oscillations are confirmed by experiments. The oscillations are permanent because 
neutrinos do not decay 

According to our approach the propagator of a particle of spin | with mass is given by formula (jTU)) and in the 
neutrino case we may deal with an f(x) as a third order polynomial thus obtaining three physical masses mi < m.2 < 
TO3, in both cases of Dirac or Majorana neutrinos 2 . The Lagrangian of interaction neutrinos/other particles can be 
divided into two parts Cf (charged current) C^ c (neutral current): 

C? c = --£= J°°W» Cf c = --J— J™Z<* (26) 

2 v2 2 cos t>w 

g electroweak coupling, 6w weak angle, W°, Z a are W^^Z®, and J^ c , Ja° charged and neutral currents respectively. 
With the neutrino masses tending to zero the interactions conserve L^^L^^Lt lepton numbers: 

L e = const L M = const L T = const (27) 

with non vanishing masses the neutrino mass term in field theory does not conserve lepton numbers. We have 

VtL = ( 28 ) 

i 

where Vi is the field of mass m, and U is the unitary mixing matrix. 

For neutrinos with definite masses there are two possibilities: if the total L = L e + L M + L T is conserved they 
are Dirac particles (four component spinors), if there are not any conserved lepton numbers they are two component 
Majorana particles 3 (no difference between neutrinos and antincutrinos). Neutrino-less double /3-decay (A, Z) — > 
{A, Z + Z) + er + e~ is forbidden if massive neutrinos are Dirac particles. 

If neutrinos are Dirac particles the mixing parameters are three rotation angles ^12,^23,^13 plus a phase factor 5; 
if they are Majorana's there are two more phases a and /3 that are irrelevant for the oscillations and play a role only 
in the neutrino-less, double-/? decay Q. More specifically we have 

C13 S 13 e~ u 

1 I • I -S12 Cia I • I e ia I (29) 



1 











C23 


523 





— 5*23 


C*23 



-W 5 Q 



23 



2 This alternative has not been established yet on experimental grounds. 

3 Processes in which the total L is conserved like (1 -> e + 7 or similar ones are permitted in case of mixing Dirac massive neutrinos 
(probabilities much smaller than the experimental upper bounds). 
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where Sij = sin^ij and Cij = cos 9ij. According to the previous arguments one easily obtains the probabilities, for 
neutrinos, of changing or not changing flavors in vacuum (see the case of atmospheric neutrinos). 

Another relevant approach deals with neutrinos in a medium (example: solar neutrinos). If we reduce our analysis 
to the transition v,, <=^ v T in vacuum we obtain 



P(i/„ -+ Vr ,L) = |K,i|^0)| 2 = C7 2 3 S 2 3 sin 2 ^-^L * C^ain 2 ^ (30) 



with Am 2 3 = m\ — m|, and Piy^ — > v^, L) = 1 — P(v^ — > v T) L). 

While the oscillations of neutrinos in the vacuum are a kincmatical phenomenon, in highly dense media it becomes 
a dynamical phenomenon because of the interaction neutrino-matter (electrons, quarks ...) mainly due to v e -e 
interaction. The nuclear fusion reactions in the core of the Sun produce electron neutrinos in a high density medium. 
Afterwards neutrinos cross a decreasing density medium before reaching the surface of the Sun [J] . 

One then can describe this situation with a phenomenological potential V(r) = \f2GpN e (r) 1 Gf is the Fermi 
constant, N e (r) is the electron density at a distance r from the centre of the Sun. Accordingly we must deal with 
neutrino effective masses and effective mixing angles; they are different in comparison with those in the vacuum. In 
particular one obtains for the 12 angles 9^ the following formula 5 

tan2^= '-t7 12 , ( 31 ) 
5^ cos 26>i2 -A 

with Sf n = m\ — m 2 , A = 2^/2GfN e E, and E is the neutrino energy 0, @. 

Numerical values: The function f(x) introduced in this and in the previous note, in the approximated form of 
a third-degree polynomial, permits to calculate masses and other physical properties of the fundamental particles, 
namely quarks, neutrinos, gluons, charged leptons ... The calculation is particularly simple, qualitatively and quanti- 
tatively significant, for the massive neutrinos (7711,7712,^,3) from the data 0] (whose values might change) where one 
can ignore the electro-weak coupling and deal with the zero order only: 

d m = ml-m\~ 76.6 ± 3.5 {meVf (32) 

|Am 2 | = 777 3 ! -i(m 2 ,-m 2 J~ 2 380± 270 (melO 2 (33) 



Indeed we must consider 



f(x) = f(x, A , Ai, A 2 , A 3 ) ~ A + Xix + X 2 x 2 + X 3 x 3 (34) 



a universal function for all weak interactions. Consequently the zeros of x — f(x) — 1 coincide with those of the charged 
leptons and for the neutrino masses 7771,777,2,7773. More precisely, if we assume on the basis of universality that the 
ratios of the masses w.r.t. the lighter one are identical in both cases, we obtain: 



777? m 2 ml ml 



, = 1 = ^ = ^ x _ = '^ = ^ x+ = ^ = ^ (35 ) 

777J 777 g 777g mj 777g 777J 

x^m\ -m\~ 76.6 {meV) 2 m x ~ 4.235 x 1(T 5 eV m 2 ~ 6 x 1CT 3 eV m 3 ~ 1.4 x lCT 1 eV (36) 
From these numbers we get 

|Am 2 | ~ 2.15 x 1(T 2 (eV) 2 (37) 

A similar calculation becomes unrealistic for quarks, and basic bosons because of the strong coupling that is present 
in addition to weak and electromagnetic forces. 



4 All neutrinos interact with electrons and quarks by neutral currents, but only electronic neutrinos interact with electrons and quarks 
via charged currents. 

5 Notice the resonance condition {9"\ = 7r/4) at the density 
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V. CONCLUSIONS 



As we proposed, the insertion of an analytic function f(x) (for x = p 2 /m 2 ) into the relativistic energy-momentum 
formula of the forceless particle preserves the relation between a certain Levy stochastic process and quantum rela- 
tivistic mechanics as in the usual case [3, 0. Furthermore the phcnomenological reduction of f(x) as a third degree 
polynomial makes the Feynman graphs convergent at any order (due to renormalizablc theories) and creates three 
poles (in the x-variable)that may be related qualitatively with three masses, and consequently with the flavor phe- 
nomenon (see details in the Section Hence we derive two propagators for quarks (charge — ^ and +§), one for 
charged leptons and one for neutrinos. 

As far as neutrinos are concerned we obtain quantitative results that are compatible with the experimental data @ ■ 
This is possible (in neutrino's case) because they couple among themselves or with other leptons via weak interactions 
only and f(x) can be assumed universal within the weak interaction scenario. Hence we can calculate the neutrino mass 
spectrum from the zero order propagator, and ignore the expansion of the renormalized weak coupling. Impossible to 
consider only the zero order propagator for other fundamental particles (such as quarks etc..) because of their e.m. 
and strong couplings. 



Appendix A 

We recall that (in one dimension for simplification) starting from the Levy-Khintchin formula for cases of centered, 
symmetric laws (as our physical case is, v{dx) symmetric measure) we have 

P 2 u 2 f 

r](u) = h / (cosux — 1) v(dx) (Al) 

2 Jr. 

u(dx) being a Levy measure, and we obtain 

id t ^{x,t) = --r ] {d x )i 3 {x,t) = -^-dl-ilj{x,t)-- f [i)(x + y,t-iP(x,t)]v(dy) (A2) 
T 2t t J r 

For stationary solutions and plane waves 

Mx,t) = e- lEot/a (b(x) a = — Mx) = e ±lux (A3) 

c 

finally we obtain equation (TJJ 

Q 



E 4,( x ) = r){u)4>{x) (A4) 



r 



Appendix B 

Remark that from (1231) we have 



v 3 

and then 



-A 2 - A 3 ± VA 

X± = 2A^ (B1) 

A = (A 2 -A 3 ) 2 -4A 1 A 3 -4A^ + 4A 3 (B2) 



A 2 A} — 1 /T\n\ 

— — = — 1 — x+ — x + — = x + + x_ + x + x- (B3) 

A 3 A 3 
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